Evidence of Critical Behavior in the Disassembly of Nuclei with A ~36 
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A wide variety of observables indicate that maximal fluctuations in the disassembly of hot nuclei 
with A~36 occur at an excitation energy of 5.6±0.5 MeV/u and temperature of 8.3±0.5 MeV. As- 
sociated with this point of maximal fluctuations are a number of quantitative indicators of apparent 
critical behavior. The associated caloric curve does not appear to show a flattening such as that seen 
for heavier systems. This suggests that, in contrast to similar signals seen for liquid-gas transitions 
in heavier nuclei, the observed behavior in these very light nuclei is associated with a transition 
much closer to the critical point. 

PACS numbers: 25.70.Pq, 24.60.Ky, 05.70.Jk 



Most efforts to determine the critical point for the 
expected liquid gas-phase transition in finite nucleonic 
matter have focused on examinations of the temperature 
and excitation energy region where maximal fluctuations 
in the disassembly of highly excited nuclei are observed 
IH El ■ A variety of signatures have been employed 
in the identification of this region ll U a J a Jj , |5L la \A la |9j, 
OH EH and several publications |EllBElHniI have 
reported the observation of apparent critical behavior. 
Fisher Droplet Model analyses have been applied to ex- 
tract critical parameters jo, E3| which are very 
close to those observed for liquid-gas phase transitions in 
macroscopic systems Data from the EOS 01 an d 

ISiS Q collaborations have been employed to construct 
a co-existence curve for nucleonic matter . Those anal- 
yses have proceeded under the assumption that the point 
of apparent critical behavior was the true critical point of 
the system 0,0- However, some recent theoretical treat- 
ments suggest that apparent signals of critical behavior 
may be encountered well away from the actual critical 
point 0, 0] and applications of A-scaling analyses have 
suggested that the observation of scaling and power-law 
mass distributions [tR fl^ | are not sufficient to identify 
the true critical point. Recently, several of the present 
authors suggested that disassembly for heavier systems 
occurs within the co-existence region, at temperatures 
well below the critical temperatures [lflj . It was further 
suggested that the decreasing importance of Coulomb ef- 
fects makes the lightest nuclei the most favorable venue 



for investigation of the critical point. 

In this letter we report results of an extensive inves- 
tigation of nuclear disassembly in nuclei of A ~ 36 ex- 
cited to energies as high as 9 MeV/u. We find that the 
maximum fluctuations occur at an excitation energy of 
5.6±0.5 MeV and a temperature of 8.3±0.5 MeV. At this 
same point, a number of indications of apparent critical 
behavior are seen. While this does not guarantee that 
the critical point has been reached, we also find that the 
caloric curve does not exhibit a plateauing at the point 
of maximum fluctuations, in contrast to experimental re- 
sults for heavier systems 11]. These observations suggest 
that the critical point for these very light nuclei may have 
been reached. 

Using the TAMU NIMROD detector [U and beams 
from the TAMU K500 Super-conducting Cyclotron, wc 
have probed the properties of excited quasi projectile-like 
fragments (QP) produced in the reactions of 47 MeV/u 
40 Ar + 58 Ni. Earlier work on systems at energies near 
the Fermi energy have demonstrated the essential binary 
nature of such collisions, even at relatively small impact 
parameters [il|. As a result, these collisions prove to be 
very useful in preparing highly excited light nuclei [j^ ]. 

The charged particle detector array of NIMROD in- 
cludes 166 individual Csl detectors arranged in 12 rings 
in polar angles from ~ 3° to ~ 170°. In these experi- 
ments each forward ring included two Si-Si-Csl telescopes 
and three Si-Csl telescopes to identify intermediate mass 
fragments (IMF). The NIMROD neutron ball, which sur- 
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rounds the charged particle array, was used to determine 
the neutron multiplicities for selected events. The corre- 
lation of the charged particle multiplicity and the neutron 
multiplicity was used to select violent collisions. In this 
work, we developed a new method to reconstruct the QP 
source. We first carry out three source (QP, Quasi- Target 
and a mid-rapidity sources) fits to the observed energy 
spectra and angular distributions of the light charged 
particles (LCP). We then employ the parameters of these 
fits to control the event — by — event assignment of indi- 
vidual LCP to one of the sources using Monte Carlo sam- 
pling techniques. We associate IMFs with Z>4 with the 
QP source if they have rapidity > 0.65y p , where y p is the 
beam rapidity. For the present analysis we have selected 
reconstructed QP events with total charge number Zqp 
> 12 from violent collisions. The 28000 reconstructed 
QP events comprise ~ 4% of violent collisions. 

The excitation energy distribution was deduced using 
the energy balance equation (2j| where the kinetic ener- 
gies of charged particles, mass excesses and average neu- 
tron contributions were considered. Using results of the 
source fits we have also evaluated (on the average) small 
corrections for undetected mass and energy (mostly as 
protons), in the sampled events. The QP source velocity 
was determined from momentum conservation of all QP 
detected particles. 

Campi plots of the natural log of the largest cluster 
charge, \nZ max , versus the natural log of the normalized 

second moment, \sxS2, [Si = ^ z '^ z """ J., ' , where nt 

is the multiplicity of QP clusters with atomic number 
Zi), are very instructive in searches for critical behav- 
ior 0, 0- In Figure 1 we present such plots for nine 
selected excitation energy bins. In the low excitation en- 
ergy region only the upper (liquid phase) branch exists. 
In the range of E* /A near 5.6 MeV/u, the liquid branch 
and the lower Z max (gas phase) branch are populated 
essentially equally. At the higher E*/A the gas phase 
branch is strongly dominant. These results indicate that 
the region of maximal fluctuations signaling a transition 
between the two phases is to be found near 5.6 MeV/u. 

To further explore this region we have investigated 
other proposed observables commonly related to fluctua- 
tions and critical behavior. Since we wish to compare our 
results with several model calculations which are based 
upon specification of the initial temperature of the ex- 
cited system, we first turn to a determination of the 
caloric curve for these A=36 nuclei. We have used two 
different techniques to derive " initial temperatures" from 
the observed apparent temperatures. The first consisted 
of fitting the kinetic energy spectra for different LCPs in 
9 different bins in E* /A with Maxwellian distributions 
to obtain the slope temperatures. To derive the initial 
temperatures it is necessary that the experimentally ob- 
served slope temperatures be corrected for effects due to 
secondary decay. We therefore employed the measured 
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(a) E* /A = 1.3-2.0 
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FIG. 1: Campi plots in different E* /A windows. 



excitation energy dependence of the multiplicity for the 
ejectile under consideration to determine initial temper- 
atures 0, E3 ■ The second technique employed results 
of a Quantum Statistical Model (QSM) calculation to 
correct the observed double isotope, H-He, ratio temper- 
atures for secondary decay effects |26j, |27| . For this work 
we employed the QSM model described in reference (2^ • 
The first technique employing observed spectral slopes 
assumes sequential evaporation of the ejectiles from a 
cooling compound nucleus source 0, |2^| while the sec- 
ond assumes simultaneous fragmentation of a reduced 
density equilibrated nucleus and subsequent secondary 
evaporation from the primary fragments |2rl 127^ . Given 
that the discussion above suggests an important transi- 
tion from liquid to gas dominance at 5.6 MeV/u excita- 
tion energy, the first method should not be appropriate 
above that energy and the second method is not appro- 
priate below that energy. The data points in Figure 2 
represent the initial temperature values determined from 
cascade corrected slopes (solid squares) and H-He isotope 
ratios (solid circles), each determined in its appropriate 
region of applicability. The two techniques lead to rea- 
sonable agreement in the transition region. We note that 
the caloric curve, defined in this manner, exhibits no ob- 
vious plateau. The temperature at 5.6 MeV excitation is 
8.3 ± 0.5 MeV. 

If the vapor phase may be characterized as an ideal 
gas of clusters then, at and above T — 8.3 MeV, 
this should be signaled by a kinetic temperature, Tki n = 
§E^ n , where Ej^ n is the Coulomb corrected average ki- 
netic energy of primary fragments. Secondary decay ef- 
fects make it difficult to test this expectation. However, 
in an inspection of the average kinetic energies for the dif- 
ferent species observed, we find that, for each E* j A win- 
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FIG. 2: Caloric curve. See text for details. 



dow, the average kinetic energy of 3 He isotropically emit- 
ted in the projectile like frame, is higher than those of 
other species. This together with simple model estimates 
indicates that the 3 He spectra are the least affected by 
secondary decay. Kinetic temperatures for 3 He, defined 
as §(-Efc — B c ) where E/, is the average kinetic energy and 
B c is the Coulomb energy, are plotted as open squares in 
Figure 2. Above T — 8.3 MeV the kinetic temperatures 
show a similar trend to the chemical temperatures but 
are approximately 1.5 MeV lower. While not perfect this 
approximate agreement provides additional evidence for 
disassembly of an equilibrated system. 

Returning to the characterization of the apparent tran- 
sition at E*/A = 5.6 MeV, we plot in Figure 3(a) 
the effective Fisher- law parameter r e / 1 extracteded from 
power law fits for Z = 2 - 7 of the QP charge distributions 
in the different E*/A windows. The results are plotted 
as a function of T/Tq, where T is 8.3 ± 0.5 MeV tem- 
perature derived for E*/A = 5.6 MeV. As seen in Figure 
1, at low E* I A a large residue always remains, i.e. the 
nucleus is basically in the liquid phase accompanied by 
some light particles. When T/Tq ~ 1, the charge dis- 
tribution shows a near power-law distribution with r e // 
~ 2.3. This value is close to the critical exponent of 
the liquid gas phase transition As T continues to 

increase, the charge distribution becomes steeper which 
indicates that the system tends to vaporize. The ob- 
served minimum in r e / / is rather broad. For comparison 
to the experimental data we also display, in Figure 3(a), 
the variation of r e / / effective predicted by three different 
models: standard statistical sequential decay as modeled 
by the code GEMINI (2^, the isospin dependent Lattice 
Gas Model (LGM) of Das Gupta et al. (using A = 36 
and Z = 16 particles in a cubic lattice with 64 sites ) 
d Hi| and Classical Molecular Dynamics (CMD) with 
Coulomb forces (2^|. The last, without Coulomb, has 
been shown to give results essentially the same as the 
LGM which is known to have a liquid gas phase transi- 



tion. Thus the addition of Coulomb interactions to that 
model is argued to be equivalent to the more difficult task 
of inclusion of Coulomb interactions in the LGM. In or- 
der to compare all parameters on a T/Tq scale, T for the 
GEMINI calculation is taken as 8.3 MeV. For the LGM 
and CMD calculations T is taken as the associated phase 
transition temperature for A = 36 and Z = 16. For the 
potential employed Ref.[Ii| these are, respectively, T = 
5.0 MeV and T = 4.5 MeV. Scaled in this way, the LGM 
and CMD results exhibit a variation of r e // vs T/Tq that 
is very similar to that seen in the experiment while the 
GEMINI calculation shows dramatically different behav- 
ior. 

In a LGM investigation of scaling and apparent crit- 
ical behavior, Gulminelli et al. have pointed out that, 
in finite systems, the size distribution of the maximum 
cluster, i.e. the liquid, might overlap with the gas clus- 
ter distribution in such a manner as to mimic the critical 
power law behavior with r e // ~ 2.2 15] . They further 
note, however, that at that point the scaling laws are sat- 
isfied 0| . We find that removal of the heaviest cluster 
from the distribution leads to distributions which are ex- 
ponential in nature over the entire energy range sampled. 
In Figure 3(b) the resultant exponential slope parame- 
ters, X e ffi are also plotted against excitation energy. A 
minimum is seen in the same region where r e f f exhibts a 
minimum. Once again the LGM and CMD models show 
similar behavior while the GEMINI calculation leads to 
very different predictions. 

For a further comparison we present in Figure 3(c) - 
3(e) experimental and calculated results for (c) the mean 
normalized second moment, (£2), from Figure 1, (d) the 
normalized variance in Z max / 'Zqp distribution, NVZ = 

(z""/^^) 01' ( e ) (^2mai) - the average atomic number 
of the second largest fragment F° r these additional 
parameters, which are often used to characterize a region 
of maximum fluctuations, the data show maxima at T/Tq 
= 1 . The LGM and CMD calculations show very similar 
behavior with some difference in absolute value. (Note 
suppressed zero points on y-axes of the plots in Figure 
3.) The GEMINI results are very different in each case. 

To recap, in our measurements for the disassembly of 
a small nucleus with A~ 36 maximal fluctuations are 
observed at 5.6 ± 0.5 MeV/u excitation energy and T 
= 8.3 ± 0.5 MeV. The fragment topological structures 
suggest the onset of a phase change. Comparisons with 
results of LGM and CMD (with Coulomb) calculations 
suggest that the critical point may have been reached. As 
pointed out in the introduction, recent theoretical treat- 
ments suggest that apparent signals of critical behavior 
may be encountered well away from the actual critical 
point [13, 0, 0, 0] , and therefore most of the parame- 
ters investigated, in Figures 1 and 3 may not be sufficient 
by themselves to identify the true critical point of the sys- 
tem. What differentiates the present work from previous 
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FIG. 3: The effective Fisher-iaw parameter (r e ff) (a), the effec- 
tive exponential law parameter (\ e ff) (b), {S2) (c), NVZ fluctua- 
tion (d), the mean charge number of the second largest fragment 
{^2max) (e). Solid squares with the error bars are experimental 
data and the lines are different model calculations as illustrated in 
right bottom corner. See details in text. 



identifications of points of critical behavior in nuclei, in 
addition to the fact that these are the lightest nuclei for 
which a detailed experimental analysis has been made, is 
the comportment of the caloric curve. For heavier sys- 
tems the points identified as the critical points appear to 
occur at excitation energies very close to those at which 
the onset of significant flattening in the caloric curves are 
observed • The reason for this flattening is still under 
discussion. It may reflect expansion and/or spinodal de- 
composition inside the coexistence region 
contrast no quasi plateau region in the caloric curve is ap- 
parent in the present study. Further, above the point of 
maximal fluctuations, the kinetic temperatures of 3 He 
ejectiles, thought to represent early emission with lit- 
tle contamination from secondary decay, show a similar 
trend to that of the chemical temperatures. Also, we be- 
lieve it is worth noting that the significance of the event 
topology in the T = 8.3 MeV region is also indicated in a 
Zipf 's law test (not shown) as proposed in [34| . Moreover, 
the analysis of A-scaling of Z max distribution, the fluctu- 
ation of total kinetic energy of QP and critical exponents 
also support the phase change around E*/A ~ 5.6 MeV 
[35| . Taken together, this body of evidence suggests a 
phase change in an equilibrated system at, or extremely 
close to, the critical point of this mesoscopic system. 
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